MZ-TH/12-17 



The Fermion Mass Hierarchy in Models with Warped Extra Dimensions and 

a Bulk Higgs. 



Paul R. ArcheiQ 



Institute for Physics (WA THEP), Johannes Gutenberg-Universitdt, 
D-55099 Mainz, Germany 



Abstract 

The phenomenological implications of allowing the Higgs to propagate in both AdSj and a class of 
asymptotically AdS spaces are considered. Without tuning, the vacuum expectation value (VEV) of the 
Higgs is peaked towards the IR tip of the space and hence such a scenario still offers a potential resolution 
to the gauge-hierarchy problem. When the exponent of the Higgs VEV is approximately two and one 
assumes order one Yukawa couplings, then the fermion Dirac mass term is found to range from ~ 10"^ 
eV to ~ 200 GeV in approximate agreement with the observed fermion masses. However, this result 
is sensitive to the exponent of the Higgs VEV, which is a free parameter This paper offers a number 
of phenomenological and theoretical motivations for considering an exponent of two to be the optimal 
value. In particular, the exponent is bounded from below by the Breitenlohner-Freedman bound and the 
requirement that the dual theory resolves the gauge hierarchy problem. While, in the model considered, 
if the exponent is too large, electroweak symmetry may not be broken. In addition, the holographic 
method is used to demonstrate, in generality, that the flatter the Higgs VEV, the smaller the contribution 
to the electroweak T parameter. In addition, the constraints from a large class of gauge mediated and 
scalar mediated flavour changing neutral currents, will be at minimal values for flatter Higgs VEVs. 
Some initial steps are taken to investigate the physical scalar degrees of freedom that arise from a mixing 
between the W5 /Z5 components and the Higgs components. 
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1 Introduction. 



In high energy physics today there exists two striking hierarchies. Firstly there exists a hierarchy between 
the two known dimensionful parameters, the Planck scale (~ 10'^ GeV) and the electroweak (EW) scale 
(~ 200 GeV). Secondly there exists a hierarchy in the observed fermion masses ranging from the top mass 
(~ 170 GeV) to the lightest neutrino mass (~ 10~^ - 10~^ eV). Although there is considerable uncertainty, 
with regard to the lightest neutrino mass, it is curious that these two apparently unrelated hierarchies should 
extend over a similar number of orders of magnitude. It has been known for some time that warped extra 
dimensions, in the context of the Randall and Sundrum model (RS) |T|, offer a compelling explanation of 
such hierarchies. In particular, the fundamental scale of the model is taken to be the Planck scale and a 
small 4D effective EW scale is generated, via gravitational redshifting, by localising the Higgs in the IR tip 
of the space. Further still, by allowing the standard model (SM) particles to propagate in the bulk, fermions 
with zero modes peaked towards the UV tip of the space, will gain exponentially small masses |l2l|3jlll. 
While this description of flavour naturally gives rise to a large hierarchy in the fermion masses, it offers no 
indication of the size of that hierarchy. 

As first pointed out in (5), this changes when one allows the Higgs to propagate in the bulk. The model still 
offers a potential resolution to the gauge hierarchy problem since it is found that in AdSs, without tuning, 
the Higgs profile is not flat, but peaked towards the IR tip of the space and would gain a vacuum expectation 
value (VEV), h{r) ~ g^ir+o-/:;-^ where k is the curvature scale of the space. Assuming order one Yukawa 
couplings, this spatial dependence of the Higgs VEV gives rise to a 'maximum' and a 'minimum' fermion 
zero mode mass corresponding to IR and UV localised fermions. So the SM fermion masses would now 
extend down, from the EW scale, by a factor of Q ' Where Q. ~ 10^^ is the warp factor related to the 
difference between the EW scale and the Planck scale. For example, if a w 0, then the fermion Dirac masses 
would range from ~ 200 GeV to ~ 10~^ eV. However the exponent of the Higgs VEV, a, is essentially a 
free parameter. Here we argue that if it is possible to find a plausible reason for why a should be close to 
zero, then models with a bulk Higgs would help to offer some insight as to why the observed fermion mass 
hierarchy extends over the range that it does. 

From a purely 5D perspective of the RS model, there is no strong argument for why the Higgs should 
have a special status as the only brane localised particle. Since, as far as the author is aware, there are 
no symmetries forbidding such bulk Higgs terms. On the other hand, in a conjectured dual theory, an IR 
brane localised Higgs would be dual to a purely composite Higgs, whereas a bulk Higgs would be dual to 
a partially composite Higgs (i.e an elementary Higgs mixing with a composite sector). The corresponding 
Higgs operator would have a scaling dimensions of2 + a |l6ir7ll8l|9l [TOl [TTl. Hence models with a bulk 
Higgs offer a convenient framework that allow for the scaling dimension of the Higgs operator to be easily 
varied. 

So the bulk of this paper will be concerned with investigating the phenomenological implications of chang- 
ing the exponent of the Higgs VEV, a. However there are a couple of related secondary questions. Firstly, 
by allowing the Higgs (a complex doublet under SU(2)) to propagate in the bulk, one gains additional scalar 
degrees of freedom that arise as a mixture of the fifth component of the W and Z fields and the additional 
components of the Higgs. Such scalar fields offer an important prediction that can help with the verifica- 
tion or exclusion of 5D models with a bulk Higgs. Hence in section [5] we take some initial steps towards 
investigating the imphcations of such pseudo-scalars. 
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A second question is related to the geometry of the space. There has been a considerable amount of work 
done on the phenomenology of the RS model, however this work has primarily focused on AdSs. At the 
very least, the RS model should include a bulk Goldberger-Wise scalar 1 12| and typically the model includes 
many more bulk fields. These bulk fields will lead to a modification in the AdSj geometry which is assumed 
to be small. An important question to ask is, to what extent is the existing phenomenological analysis robust 
against small modifications to the geometry? With this in mind, in addition to a pure AdSs geometry, we 
shall also consider a class of asymptotically AdSj spaces that have arisen from a scalar plus gravity system 
|[T3l . These spaces are of additional interest since they can, for certain regions of parameter space, result in 
a significant reduction in the constraints from EW observables l[T4l ITSl IT6l and flavour physics lITTl . 

Aside from these two additional questions, the central result of this paper is that there are a number of both 
theoretical and phenomenological motivations for considering a close to zero to be the optimal value. These 
include: 

• The four dimensional effective Higgs potential receives an additional positive |Op term proportional 
to a which can, for large values of a, dominate over the negative |Op term and result in an EW phase 
transition no longer occurring. 

• If the space is not cut off in the IR, then the Breitenlohner-Freedman bound implies a ^ ifTSll . 

• The holographic method is used to demonstrate, for generic geometries and generic Higgs potentials, 
that the flatter the Higgs VEV (i.e. the smaller the value of a) the smaller the contribution to the 
Peskin-Takeuchi T parameter [19]. 

• For the RS model, large values of a can give rise to strongly coupled pseudo-scalars and potentially 
large constraints from pseudo-scalar mediated flavour changing neutral currents (FCNC's). 

• It is found, for all spaces considered, that the modification to the W and Z profile will be minimal 
for smaller values of a. This will likely result in a reduction in many of the constraints from flavour 
physics. 

• Assuming anarchic Yukawa couplings, there is a general trend that suggests that the smaller the value 
of a the more the SM fermion profiles will be peaked towards the UV tip of the space. This will again 
result in the reduction of constraints from FCNC's lfT7ll20ll2T]| . 

• In order for a 4D Higgs operator, O, with a scaling dimension 2 + or, to offer a resolution to the gauge 
hierarchy problem then the corresponding 0^0 operator must not be relevant, which implies a ^ 

m. 

The outline of this paper is as follows. In section|2]the geometry and the Higgs sector is outlined. In section 
[3] the fermion mass hierarchy is discussed. In section [4| the EW constraints are computed. In section [5] 
the pseudo-scalais are investigated. In section [6] the implications for flavour physics are discussed and we 
conclude in section |7] The majority of the core equations have been derived, for a generic space, in the 
appendix. This paper should be considered as an extension of a number of existing studies on bulk Higgs 
scenarios, including Il22ll23ll24ll25ll26l . 
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2 The Model. 



This paper seeks to investigate the most minimal extension of the SM in which the Higgs propagates in the 
bulk of a 5D warped extra dimension. Hence we shall consider a bulk SU(2) x U(l) gauge symmetry and 
not the extended custodial symmetry fT/]. The fermion content remains the same as the SM and we include 
a bulk Higgs, which is a doublet under the SU(2) symmetry. As in the RS model, we compactify the 
space over a /Z2 orbifold in which the space is cut-off at the two fixed points, ruv = and rjR = R. The 
Higgs sector is then described by 

S=J d'x^[\DMO\^-Vm] + J ^/^^V^[-Vir(0)],=« + J d^x^/^[-Vuym],^o (1) 

where Vir/uv are the Higgs potentials localised on fiat branes located at the orbifold fixed points. While 
G and giR/uv are the determinant of the bulk metric and the induced brane metrics. At this point we must 
make some simplifying assumptions. Firstly we have neglected brane localised kinetic terms and secondly 
here we consider potentials of the form 

Vm = Mlm^ Vir(O) = -Mir|(D|2 + ^ir|(D|4 VuvW = Muv|0|^ (2) 

We have not included terms in the bulk or on the UV brane. In RS type scenarios, the fundamental 
coefficients of such operators are assumed to be at the Planck scale ~ 0{k^^). However, in the 4D effective 
theory, would be warped down to the Kaluza-Klein (KK) scale ~ 0{M^^ while the corresponding UV 
operator would remain at the Planck scale. A bulk l^l'^ term is assumed to be suppressed by an intermediate 
scale. Although including such a term would result in the Higgs VEV being the solution of a nonlinear 
differential equation and so result in a significantly more complicated model. 



2.1 The Geometry. 

With these assumptions the model is largely completely defined and the only remaining input is the geometry 
of the space. Here we shall consider spaces of the form 

ds^ = e-^^^'^rfdx^dxy - d? (3) 

where = diag(l, -1,-1,-1) and < r < /?. In the original RS model just gravity propagated in the bulk 
and the solution considered was a slice of AdSs, 

A(r) - kr. (4) 

However it was quickly realised that in order to stabilise the space one must also include an additional bulk 
Goldberger-Wise scalar |[T2ll . Also, in order to generate the fermion mass hierarchy and suppress EW and 
flavour constraints, the model was extended to allow the SM particles to propagate in the bulk EISIH. The 
back reaction of such bulk fields will typically lead to a deviation, from the AdS geometry, in the IR tip of 
the space. Although the space should be asymptotically AdS in the UV|^ While it is not certain how large 

^The 5D generalisation of BirkhofF's theorem i28il ensures tiiat solutions with no (or constant) bulk fields must be AdS. In the 
RS model the only fields peaked towards the UV are the light fermion modes which are assumed to have a negligible coupling to 
gravity. The other fields are typically vanishing in the UV. Hence it is anticipated that the space should be asymptotically AdS in 
the UV. 
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this IR deformation should be, it is important to ask how sensitive any result, based on AdSs, is to such 
modifications in the geometry. So in this paper, we shall also consider a class of geometries that have arisen 
from solving a scalar plus gravity model, giving solutions of the form lfT3l[T5]| . 

A(r)....lln(l-^). (5) 

Firstly, it should be noted that scalar gravity systems generically give rise to a singularity [29], here located 
at r = R + A. While it is possible to construct models that incorporate this singularity, such as soft-wall 
models, here we shall impose an IR cut-off, on the space, before it reaches the singularity (i.e. let r e [0, R] 
and require A > 0). The parameters v and A are determined by the scalar potential and boundary conditions, 
although here we shall treat them as free parameters. Note AdSs is regained by sending v and A to infinity. 

At this point it is useful to define the warp factor and KK scale to be 



o A(R) A i.^ d5A{r)\r=R 

Q = ' and Mkk = 



KK ^ (6) 



In order for the space to offer the potential for a non-supersymmetric resolution to the gauge hierarchy 
problem, it is required that the space can be stabilised, such that Q ~ lO'^, and also that Mkk ~ (9(TeV) is 
phenomenologically viable. 



2.2 The Higgs VEV. 



In the above model, a bulk Higgs would gain a VEV, (O) - ^ j' where h{r) satisfies 

dlKr) - 4 55A(r) dih{r) - M|/i = 0, (7) 

with the consistent boundary conditions being either /i|,=o,R = or 

[55/j - Muv/j].=o = and [d^h - Mi^h + Ai^h\^^ = 0. (8) 

Since we are not considering an EW symmetry breaking bulk potential we must choose the latter 'non- 
Dirichlet' boundary conditions. In other words here EW symmetry is broken on the IR brane. In the early 
days of the RS model, it was argued that models with a bulk Higgs required fine tuning in order to achieve the 
correct W and Z masses ||30l . However this work assumed a Higgs VEV that was constant with respect to r. 
From the above relations, it can be seen that a constant Higgs VEV can only be achieved by either breaking 
EW symmetry on both the branes and the bulk, such that d(i)V{<^)\i^=h = d(^Vui{(^)\(^=h = -d,^Vi]Yi^)\o=h - 
0, or alternatively breaking EW symmetry just in the bulk and forbidding the existence of brane potentials. 
Clearly the first option is finely tuned and it is not clear how the second option could be achieved. Hence 
here we would suggest that a Higgs VEV peaked towards the IR is a more natural scenario, for which the 
arguments of 1,30,1 are not applicable. 

For the RS model Q, it is then straight forward to solve for the Higgs VEV (see for example Il6ll23]|) 

h{r) = Nhe^^' (e"'"' + Be-"'"') (9) 



4 



where B and Nh are constants of integration fixed by tlie boundary conditionqj 

2k + ak- Muv , 2 _ (2^^ + Q-^ - Mir)Q2+" + B{2k - ak - MuOQ- 

^ ~ ~2k-ak-Mvy ~ ^iR(Q2+o- + BQ2-a)3 

We have also introduced the Higgs exponent, which will prove to be an important parameter, 



(10) 



~ k 



(11) 



Clearly, one can always take the limit or ^ 00 in order to gain an IR brane localised Higgs. As shall be 
demonstrated in the next section, in order to explain the observed fermion mass hierarchy it is required that 
a is small, a ~ 0, which clearly requires that M| < 0. In AdS space, without an IR cut off, a negative mass 
squared term is permitted provided it satisfied the Breitenlohner-Freedman bound ifTSl . In particular, in order 
for the total energy of the scalar to be conserved (as well as to allow for a valid global Cauchy surface) it is 
necessary for the energy-momentum flux to vanish at the AdS boundary. For AdSs, this implies M| ^ -4k^ 
and a ^ 0. Technically this bound does not apply here since we are cutting the space off before we reach the 
AdS boundary. Although the RS model is not a UV complete theory and so it is feasible that such a bound 
is applicable in a more fundamental realisation of the model. Further still, it was pointed out in Il22l that 
such a tachyonic mass term could naturally come from the Higgs coupling to gravity ~ ^'R|Op, where ^ is 
an (9(1) coupling and the Ricci scalar is = -20k^. 

While this paper will primarily focus on the 5D theory, it is worth making a brief comment about the con- 
jectured dual theory. Although we refer readers to @ [71 [8l |9l [TOl \TT] for a more comprehensive discussion. 
A bulk Higgs in the RS model is conjectured to be dual to an operator, O, of a broken conformal field theory 
mixing with an elementary source field (pQ = O|r=o- In other words it is conjectured to be dual to a partially 
composite Higgs. By computing the bulk-brane propagator and taking the large Euclidean momentum limit, 
it is found that the scaling dimension of the Higgs operator is given by 2 -1- a ||6j [TQ. Further still it is 
found that, below the cut-off A ~ k, the coupling between the source field and operator is approximately 
-C4Z) ~ ^^fpoO, where coisa. dimensionless constant. Hence the coupling is relevant for a < 1 and therefore 
a determines the level of mixing between the source field and the composite operator. The point we wish 
to emphasize here is that the parameter a determines both the scaling dimension of the Higgs operator and 
its level of compositeness, i.e. the larger a the more composite the Higgs. A consequence of this result, 
discussed in [6|, is that in order to resolve the hierarchy problem, one require that the corresponding 4D 
Higgs mass operator, O^O, is not relevant. This again implies that a ^ 0. 

The statements of the previous paragraph applies to AdS 5 and while one would anticipate that they should 
approximately hold when one modifies the geometry, due to the lack of an analytical expression for the 
bulk-brane propagator, this is challenging to verify. So for the remainder of this paper we shall focus purely 
on the 5D theory. Turning now to the modified metric Q, the Higgs VEV can again be found by solving ^ 
with ([8]) to get 



{v^-2)a + 4 2v^-4 



, 2ak{R + A-r) 




(12) 



^In the special case when or = then h(r) = Ni,(e^'"' + r(Muv - 2k)e^'"') and 

2 Muv(l +21nfi) + MiR(l -21n£l)-4ytln£l + i?MuvMiR 



n*Aif.(l+R(Mvw-2k)y 



5 



where U and M are confluent hypergeometric functions or Kummer functions and 

{{2k -ak- M^y)(R + A) + ^) t/ (^^^, ^^,2ak(R + A)) + ^-^,2ak(R + A)) 

{(2k -ak- M^y)(R + A) + ^) M 2a^(/? + A)) + ^-^-iz^ M {^-^ , ^^,2ak(R + A)) ' 

(13) 

Bearing in mind that for large z > 0, M{a,b,z) ~ e^f~^{l + <9(z"')) and U{a,b,z) ~ zr^il + 0{z'^)), then 
far from the IR tip of the space 

n; / 2(g+2) , 2(tt-2) \ 

hir) ~ A^/^f-^*^'" ( (/? + A - r)"^e"''^'' + B(R + A- rf^e-"'"'] . 

In other words, in the UV tip of the space one just finds a small power law correction to (|9]). However, this 
approximation breaks down in the IR, particularly for v < 1 when "^,2"'"'"^ is large. This is a potentially 
quite interesting region of parameter space in which the Higgs VEV typically grows more than exponentially 
towards the IR. Unfortunately this greater than exponential growth makes carrying out the numerical studies, 
conducted later in this paper, quite challenging. So we leave a thorough study of this region of parameter 
space to future work. Even when v > 1 , the above approximation is still not very robust in the IR. In practice 
we find a better approximation is 



hir) « Nhe^^''"*'' (R + A- r)^^ K,^ 

j,2 

where is a modified Bessel function. Note in practice it is found that B is negligibly small for much of 
the parameter space. This raises the important point that, since we are breaking EW symmetry in the IR and 
the Higgs VEV is peaked towards the IR, the majority of results of this paper are relatively insensitive to the 
value of Muv- This further supports one of our initial assumptions, notably the exclusion of a lOI'* term on 
the UV brane. For the remainder of this paper we shall fix Muv - k. 



2k{{v^ - 2)a + 4){R + A - r) 



(14) 



2.3 Fixing Some of the Input Parameters. 



Even with this relative insensitivity to the UV potential, one may be concerned with the number of free 
parameters introduced. In particular, in addition to the warp factor and KK scale, the modified metrics have 
the additional geometrical input parameters v (the IR modification to the curvature) and A (the position of 
the singularity). While a bulk Higgs, with the assumed potentials (|2]), gives rise to an additional four free 
parameters. Two of these parameters are determined by fitting to the EW scale and the Higgs mass. This 
leaves four free parameters in addition to those of the RS model with a brane localised Higgs. Here we 
would argue that this enlargement of the parameter space has arisen from relaxing assumptions of the RS 
model that (from a bottom up perspective) are questionable. 



We will now move on to fix the two parameters Mir and Aj^. Throughout this paper we shall make the 
working assumption that the Higgs Mass is 125 GeV In practice, provided that the Higgs is lighter than the 
KK scale, changing the Higgs mass will not significantly change the results. The SM Higgs is taken to be 



the lowest mass eigenstate of (48 1 with the boundary conditions 



and 



dsH - Mmih + H) + Air -h^H + h 



0. 



(15) 



r=R 



where h and H are defined in (46 1. Of course one could also impose Dirichlet boundary conditions (DBC's) 
on the Higgs. This would significantly change the phenomenological imphcations of this model and would, 
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for example, result in no light Higgs. However in light of recent LHC results we consider this scenario 
disfavoured. Also note that if the Higgs was just localised on the IR brane then the first 'c^j' terms in ^ 

and ( 15 1 would not be present. This would imply the familiar relation h 



Higgs mass being given by just Mir. 



{iS^)' ^^'^ result in the tree level 



For the RS model, ( 48 1 can be solve analytically to get 

/r = Nne^"- {ja I -^^^ \+/3Yj '^^^ 11 (16) 



(H) kr\ ( (H) kr 

m„ e I I m„ e 



where we have made the usual KK decomposition, H{x,r) = Y^nfn^\^)tl^"'Kx), such that d/jd^H^"^ = 
-m„ ^H^"\ While J and Y are Bessel functions and /? is a constant of integration fixed by the boundary 
conditions. When considering the modified metrics one must work with a numerical solution for the Higgs 
profile. In addition to the Higgs mass we also fit to the Fermi constant, Gf - 1.166367(5) x 10"^, the Z 
mass,Mz ^ 9 1.1 876 ±0.002 1 GeV and the fine structure constant, q'(Mz)"^ = 127.916 ±0.0015 EU. These 
are given, at tree level, by 



2 _j_ fyf 2 



■(r) 



where fl^'^^ are the photon and W profiles, defined in the appendix. While yj^^'^'''-^^''''''-^ ^re the fermion 



zero mode profiles for fermions with a bulk mass parameter c^, given in (70 1. For the purpose of this fit, we 
assume that the muon, electron and neutrinoes all have a universal bulk position cl - -cr = 0.7. Also note 
that, with the absence of a right handed neutrino zero mode, the charged pseudo-scalars will not contribute 
at tree level to the Fermi constant (see section |5). 



The results have been plotted in figure [T] This brings us to the first motivation for favouring small values 
of a. In the standard Higgs mechanism one requires a negative |<l)p term, in the Higgs potential, in order to 
generate a non zero Higgs VEV and break EW symmetry. When one allows the Higgs to propagate in the 
bulk of an extra dimension, the 4D effective potential will gain, at tree level, a positive contribution to the 
lOp term loosely related to d5<^. The larger the value of a, the larger this contribution will be and if a is too 
large then EW symmetry will no longer be broken at tree level. This can be seen explicitly in figure [T] As 
a is increased, the require negative Mir also increases, in order to compensate the positive contribution and 
at some point {a ~ 5 for the RS model) the positive contribution dominates and EW symmetry is no longer 
broken. 



One should be cautious about taking the values of a in figure [T| for which this happens, too literally. Where 
it is not possible to find values of Mir and Air that break EW symmetry, one should of course consider 
additional operators, in particular a bulk l^l^^ term and it is not clear to what extent this effect continues 
beyond tree level. Rather figure [T] gives an indicator of the values of a for which breaking EW symmetry, 
in a fashion that is compatible with EW observables, becomes difficult. What is surprising is that for certain 
geometries, in particular geometries with small v and large A, this can happen for relatively small values of 
a. 
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(a) The RS Model 
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(b) Modified metric with v = 10 



(c) Modified metric with v = 3 



Figure 1: The values of Mir and Air after fitting to Mz, O!, Gf and - 125 GeV. For the modified 
metrics we have taken = 0.5 (dash-dot), kA = I (solid), kA = 1.5 (dash-dash) and kA = 2 (dots). We 
have fixed O = 10^^ Muv = k and Mkk = 2 TeV. 

3 Fermion Masses. 



Having found the form of the Higgs VEV, we can now move on to look at the fermion masses. The fermions 



will gain a mass via the Yukawa coupling, Y, to the Higgs. I.e. in addition to the terms in (68 1 the action 
will also include 

where *F is a doublet under SU(2) and X is a singlet. If we split the spinor into its chiral components 
X = XL + Xr arid likewise for ^ (as in the appendix) then, after compactifying over the S ^ /Z2 and choosing 
appropriate boundary conditions, only ifri and xr will gain zero modes. So a low energy chiral theory is 
achieved with T s {L, Q] and X s {e, u,d}. It can now be seen that, due to this mixing by the bulk Higgs, 



the profiles given in (69 1 will not be mass eigenstates. To compute the mass eigenstates, if we define 



An,m) _ _j_y r 

" V2 J 



drh/r'f: 



and 



Mn,m) _ _j_y r 



drhfr'fr. 



(18) 
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then the fermion mass matrix will be given by 






( y(0,0) 
y(l,0) 



y(2,0) 
^LR 





y(l,l) 
^RL 





y(0,l) 
^LR 
y(l,l) 
^LR 

y(2,l) 
^ LR 







y(2,l) 
^RL 

w 



( y(0) \ 

4;; 



(19) 



Note we have restricted this discussion to one generation, in reality Ylr and Yrl will be block 3x3 matrices. 
This matrix can then be perturbatively diagonalised (see for example ||32| ) such that the zero mode mass 
eigenvalue is found to be approximately 

y(0,M)y(n,m) y(m,0) 

n,m=l v'"n )\"^m ^ i^r ) 

While the profile of the mass eigenstate of the^^ zero mode, for example, will be 

y(0,«) 

h Zj. Or) y(0,0)/« ^ 1.^^^ 

n=l V^n - J^LR ) 

So the SM fermion masses will be given by y''^'^^ with an 0{m^/M^^) corrections and will have the profiles 

/q^'^^ with (3(mo/MKK) corrections. In this paper we will be largely concerned with the lighter fermions and 
so we will neglect these corrections. However these corrections should not be completely forgotten since 
they may lead to phenomenological deviations from the brane Higgs scenario, particularly with regard to 
top physics. It is also possible that in models with extended Higgs sectors, such as models with a custodial 
symmetry [27] or models with axial gluons [33.1 . these corrections may be enhanced. 



For the RS model it is then straightforward to analytically integrate over ^ and the fermion profiles ( 70 1 
giving 



(1 - 2cl)(1 + 2c r) ((Mir - 2 - a)a'' + B(Mir - 2 + a)a-<') 

— Mkk a 

V2 } (Qi-2^'. - l)(ni+2'^« - 1) Air (Q'' + BD.-")^ 



+ B , (22) 



^LR - — ^1 



i2-CL + CR + a) (2-CL + CR-a) y 

where we have introduced the 'assumed (9(1)' coefficients Mir,uv = Mjr uv^ \ -^ir - ^ir^^ and Y - Y ^/k. 
Although in practice, since the EW scale is not the same as the KK scale, these parameters are not (9(1) (see 
figure [T]l. This should be contrasted with the analogous expression for the brane localised Higgs (see for 
example [[341), 

y(0,0) _ Y j (1 - 2cl)(1 + 2cr) ^j-c,,+c, .... 

""^ '\(al-2^^-l)(Ql+2c._l)" (23) 

where the 4D Higgs VEV is V4 ~ 246 Ge\n The range of fermion masses have been plotted in figure pi 
For the modified metrics, Y^ does not have a neat analytical form but the resulting distribution of fermion 



masses is similar to that which is plotted in figure [2(b)l 



''In practice V4, in the RS model with a brane locahsed Higgs, will still differ from the SM value. See t35J . 
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(a) Brane localised Higgs 



(b) Bulk Higgs {a = 0.01) 



y(0,0) 

Figure 2: The approximate range of fermion zero mode masses, -y-, for a bulk and brane Higgs in the RS 
model. The fermions are peaked towards the UV (IR) when cl > 0.5 ( < 0.5) and cr < -0.5 (cr > -0.5). 
n - 10^5 and Mkk = 2 TeV. 



This brings us to the central result of this paper and the initial motivation for this study. It is already widely 
known that warped extra dimensions offer a potential explanation for the large hierarchies that exist in the 



observed fermion masses [3l IH. As can be seen in figure 2(a) order one changes in the bulk mass 
parameter result in exponential changes in the zero mode mass. However the model offers no indication 
of what the overall size of the mass hierarchy should be. When we consider a bulk Higgs this situation 
changes. When the fermion profiles are peaked towards the IR brane, the Qi-CL+c«+a tej-m will dominate 

and niQ ~ FMkk^jr ~ 200 - 300 GeV. Although when the fermions are peaked towards the UV brane the 

second term will dominate and mo ~ FMkk^jr f^"^ "^- Hence if we assume there are no hierarchies 
in Y and we do not have a split fermion scenario, then the difference between the lightest and the heaviest 
fermion Dirac mass term is < 

These are of course approximate relations, but the calculated values are plotted in figure [3] This figure re- 
quires a number of remarks. Firstly, what has been plotted is not the absolute maximum and minimum mass 



but rather the location of the upper and lower 'plateaus' in figure 2(b) These plateaus are not strictly flat 
but increase (or decrease) logarithmically and also the bulk mass parameters, for which the lower 'plateau' 
begins, increases as a is increased (see analogous plots in 1211). We have also assumed F w 1. By naive 
dimensional analysis it is suspected that perturbative control of the theory is lost if 7 > 3 1361, but there 
is nothing, except naturalness arguments, forbidding Y from being small. Likewise we have rejected the 
possibility of a split fermion scenario in which the left handed fermions are peaked towards the UV and 
the right handed fermions are peaked towards the IR or vice-versa. Such a scenario would undoubtably 
give rise to large constraints from FCNC's and so are typically disfavoured by fits including relevant flavour 
observables (see section[6]). 
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1 o-3°l > > > > > > > > ^ 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

a 

(a) 'Minimum mass' =7^'''/? with = 5, = -5 




(b) 'Maximum mass' = yj^''V? with Cl = -3,cr = 3 



Figure 3: The approximate upper and lower limits of the fermion zero mode mass, assuming Y » I and 
one does not have a split fermion scenario. Here we have considered the RS model (black) as well as the 
modified metrics with v = 10 (blue), v = 5 (green) and v = 3 (red) with = 1.5 (dash-dash line), kA = I 
(solid line) and kA = 0.5 (dash-dot line). We have set Muv - k, Mkk - 2 TeV and Q. = lO'^. The remaining 



input parameters have been fixed as described in section 2.3 The 'maximum fermion mass' does decrease 
for large a and it is suspected that it should flow towards ~ 174 GeV as a ^ oo for large Mkk- These large 
values of a have not been included in the plot since they would not give rise to a phenomenologically viable 
electroweak phases transition, as described in the previous section. For large values of a the lower 'plateau' 
will not have begun for the bulk mass parameters cl = ^ and cr = -3 and so the minimum fermion mass is 
evaluated at cl = 5 and cr = -5. 
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If we accept these assumptions, then we arrive at a potentially interesting result, first pointed out in O. 
For spaces close to AdSj, with a small value of a < 0.1, the lower bounds on the fermions masses are 
remarkably close to the suspected value for the lightest neutrino mass, 10~^ - 10~^ eV|^ This alone does not 
explain the size of the neutrino masses, since the above results are for the Dirac mass term. Hence one must 
also demonstrate why either the Majorana mass term is forbidden, or alternatively explain why it exists at 
approximately the same scale. There are many possibilities for achieving this, for example, the latter option 
could be achieved if the origin of the Majorana mass term had a similar profile to the Higgs VEV. However 
these possibilities would necessarily require an extension of the minimal set up considered in this paper and 
so here we leave them to future work. Also, a bulk Higgs scenario offers no explanation as to why there 
is a six orders of magnitude 'desert' between the electron mass, and the heaviest neutrino mass or why the 
leptons are lighter than the quarks. In other words, this work is concerned with the range of the Dirac mass 
term, while statements concerning neutrino masses necessarily require an extended model. 

Having said that, here we would argue that if one could understand why a should be small, then models 
with a bulk Higgs would offer some insight into why the fermion mass hierarchy is the size that it is. Having 
already offered two motivations for favouring a small value of a, we shall now move on to demonstrate that 
many of the phenomenological constraints on warped extra dimensions will be at their minimum for smaller 
values of a. 



4 Electroweak Constraints. 



One of the first tests of any model, that modifies the Higgs mechanism of the SM, is that of being able 
to suppress corrections to the precision EW observables. Typically it is found that such corrections are 
dominated by modifications to the W and Z propagators and are often parameterised in terms of the 'oblique' 
Peskin-Takeuchi parameters llT9l . Such a parameterisation assumes that the corrections to gauge-fermion 
couplings are negligible and the new physics can be described by the effective Lagrangian 

Xeff . = - [p^ -ml- Uzip)) - IV; (p2 - m\ - Xl^{p)) + . . . (24) 

where p = pi^ is the four momentum. The effective Lagrangian also contains terms for the photon as well as 
terms mixing the photon and the Z, but these corrections are found to be zero, at tree level, in models with 
extra dimensions. Expanding the H's as 

uip) = n(0) + p^u'iO) + ^n"(0) + . . . 

then the Peskin-Takeuchi 5 and T parameters are defined to be 

aT = luziO) - — nw(0) I and aS = 4 sin^ e„ cos^ 6^ (u'^iO) - n;(0)) , (25) 
in^' \ sin 0m j 

where 6„ is the weak mixing angle. It is possible to introduce additional parameters in order to account for 
the higher order corrections [38]. However here it is found that, with the absence of a custodial symmetry 
Il27il . the constraint is dominated by contributions to the T parameter 

^While there is considerable uncertainty in this mass, assuming a normal hierarchy, the recent measurement of a large value of 
9\i would probably favour a lighter minimum neutrino mass. See for example 1371 
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In order to compute the corrections to the propagator here we shall largely follow ifTSl and expand the 5D 
field in terms of the holographic basis. This method is completely equivalent to a KK expansion and does 
not rely on the existence of a dual theory. In the interests of generality, we shall also work with the generic 
metric considered in the appendix (43 1. In particular we factorise the 5D fields 

Z^{p, r) = G^^\p, r%{p) such that Zf,{p, ruv) = Z^{p) 



(26) 



and likewise for the W. We also impose Neumann boundary conditions (NBC's) on just the IR brane, 
d^Zf^li^ = 0. Bearing in mind that G^^'> will satisfy dsia'^b-'^dsG^^^) - o^bMlO^^^ + bp'^G^^^ = 0, the tree 
level effective Lagrangian will be given purely by the UV boundary term 



Xeff. = ^Z, [G'^^a^b-'d5G^%_^^ Z^ + W; [G^^Vb-'dsG'"^^] 



W1 



Jr=ruv 



(27) 



We have once again neglected possible brane localised kinetic terms which are known to reduce the EW 
constraints [,39,1 . To proceed further we define 



Pw,ziP^ r) 



a^b-^d5G^^'^\p,r) 
G(^-^\p,r) 



(28) 



such that 



d5Pw,z + a-^bP'iy - a^Mly + bp^ = 0. 



We can now match P]v,z to the oblique parameters by expanding in four momentum 



(29) 



PwAP^ r) = P%\{r) + p'P'y,'y{r) + - p'' P'^:^,ir) + ... 



.2 d(1) 



-P^P^^^ V 

2' w,z'^ 



(30) 



and equating ( 24 1 to (27 1 giving 
niv,z(0) = P 



m I 

W,Z\r=r]j\ 



m 



w,z 



and U'^^^iO) = pZ\1__ + 1 



W,Z\r=rD\ 



(31) 



Where P"^ can be found by expanding (|29[) into a set of coupled equations 



a^bMl/2 = 



dsP%]^ + la-^bP^PZ, + b^0 



w,z^ w,z 



^5 W,Z ^'^^ W7^ W7 ^ W7 / 



W,Z^ w,z 



w,z- 



(32) 



As already mentioned, the greatest contribution to the EW constraints comes from the T parameter which is 
given by P^^lz\r-n]^/^ Without loss of generality, we can set = 1, <3(riR) = Q.'^ and fl:(ruv) = 1- Then using 

the IR boundary condition, PvKz(nR) 
the T parameter) as 



p'-^^(riK) - 0, (32 1 can be rearranged to give 



W,Z If— ruv 



(and hence 



3(0) 2 1 
W,Z 



r=ruv 



'-'5^w,z\ 



r=nR 



"5^w,z\ 



r=ru\ 



+ M\ 



w,zi''=njv 



^w,z\r=nK- 



(33) 



It should be stressed that this holds for generic Higgs VEV's and generic 5D geometries. This then brings us 



to the point of this discussion. Bearing in mind the definition of Miy,z (49 1, it is now straight forward to see 



that the 'flatter' the Higgs VEV (i.e. the smaller the value of a), the smaller the last two terms of (33 1 will 
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be and the smaller the contribution to the T parameter will bqj This should not be a particularly surprising 
result. One gets tree level corrections to the EW observables because the Higgs mixes the SM gauge fields 
(the zero modes) with the KK gauge fields that are peaked in the JR. The less the Higgs is peaked towards 
the IR then the weaker this mixing will be. Alternatively, in the language of a possible dual theory, the less 
composite the Higgs, the less it will mix the composite states with the elementary states. 




12 3 4 
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Figure 4: The constraints on Mkk (not the mass of the first KK gauge boson) from EW precision tests for 
the RS model (black) as well as the modified metrics with v = 10 (blue), v = 5 (green) and v = 3 (red). We 
have used kls.- I (sohd line), M - 0.5 (dash dash), mn - 125 GeV, Muv - k and fixed Q. = 10^^. At each 
iteration we have fit to a, Mz and G / and compared with the latest S - T ellipse ( 34 ) at 95% confidence level. 
For = 0.5, such an analysis is plagued by small numerical inaccuracies which results in not completely 
smooth curves. 



The contribution to the first two terms of ( 33 1 are quite sensitive to the geometry of the space. It was found 
in l[T4l[T5l[T6l . that the contribution to these two terms could be significantly reduced in the modified spaces 
(|5]l with V < 1 and kA <, I. However, as discussed in section 2.2 this region of parameter space typically 



results in the Higgs VEV growing more than exponentially which would increase the contribution from the 



last two terms of (33 1. So a minimum in the EW constraints can be found by introducing a moderate tuning 
such that the Higgs VEV still grows exponentially and v < 1. 



Before we can compute the size of these constraints, we must first perform a field redefinition in order to 
absorb the non-oblique corrections to the gauge-fermion couplings, 



W, 



g4 



*This should not be mistaken for a proof since we have not demonstrated that the a dependence of the last two terms dominates 
over the first two terms. It is probably possible to contrive a space for which this is not true. However, one can see that it will hold 



for all spaces in which P'jJ^ is growing towards the IR and d^P^!!^^ oc a" with « > 0, 



i(0) 
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RS 


V - 10 


v = 5 


V = 3 


0.5 


2.45 


2.23 


1.68 


0.87 


1 


2.45 


2.27 


1.82 


1.07 


1.5 


2.45 


2.30 


1.89 


1.18 



Table 1: The mass of the first gauge boson relative to the KK scale, m. '' /Mkk- 



and likewise for the Z field. Here g4 is the 4D effective coupling found by fitting to ( 17 ). This is only possible 
if one assumes universal fermion couplings and hence we again assume a universal fermion position with 
cl = -cr - 0.7. In practice this is not a good approximation for the heavy quarks. Although the modified 
metrics, with more realistic fermion positions, have been studied in llT6ll and it was still found that certain 
regions of the parameter space still have significantly reduced constraints. 



The constraints on the KK scale can then be calculated by comparison with a fit to the S and T parameters 
(assuming a vanishing U parameter) POl . 

S = 0.07 ± 0.09 r - 0.10 + 0.08 Pcorrelation = +0.88. (34) 

These constraints have been plotted in figure |4] Note these constraints are for the KK scale and not for mass 
of the first KK gauge field. The masses of the first KK photon (or gluon) are given in table [T] Here we 
find that relatively small shifts in the geometry can result in significant reductions in the mass eigenvalues 
relative to the defined KK scale. It is believed that this is partly responsible for the reduction in the EW 
constraints found in |[T4l[T5l[T6]| . 



5 The Pseudo-Scalars. 



Having demonstrated that EW constraints will generically be at their minimum for smaller values of a, we 
shall now move on to look at how models with a bulk Higgs can be potentially falsified. Before spontaneous 
symmetry breaking the model considered here contains four 5D massless gauge fields (4x3 transverse 
degrees of freedom) and the Higgs, a complex doublet (4 scalar degrees of freedom). After compactification 
and the breaking of EW symmetry the model contains four 4D massive gauge fields (4x2 transverse degrees 
of freedom, plus 4 longitudinal degrees of freedom) and a Higgs particle (1 scalar degree of freedom). Hence 
the model must necessarily also include an additional three scalar degrees of freedom. In this section we 
shall demonstrate that, for a ~ 0{\) or less, such pseudo-scalars will gain masses at the KK scale and hence 
the observation of a Z' or W at the LHC should be associated with the existence of a corresponding charged 
and neutral scalar. 



These pseudo-scalars have been previously investigated for warped spaces in 10311411 and are well known in 
models with universal extra dimensions. As found in the appendix, these pseudo-scalars arise as a mixture 
of the A5 component of the gauge field and the Higgs components, tt,-, see (58 1 and (63 1. Note the larger the 
value of a the more these scalars will be dominated by the Higgs components tt; and hence more of the A5 
component can be 'gauged away'. 
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5.1 Consistent Boundary Conditions. 



If we begin by considering a gauge field in the 5D space ( 43 1 then variation of the action yields the boundary 
term 

"1 
2' 



= 0. 



In order for such a boundary term to vanish one must impose either NBC's on the components and DBC's 
on the A5 components or vice versa. If one is considering just an interval then one can, in theory, impose 
non-trivial boundary conditions. However if one is compactifying the space over a S ^ /Z2 orbifold, such 
non-trivial boundary conditions are forbidden by gauge invariance across the fixed point^ Hence if one 
imposes NBC's on the A^ field, in order to gain a zero mode that is associated with the SM field, then one 
must impose DBC's on the A5 component. 



As for the Higgs particle, the consistent boundary conditions for the Higgs components are either DBC's, 



i'-=nR,njv 



^ 0, or non DBC's 



\b ^dsTTi- MiRjTi + AiRh^jTi] =0 and \b ^dsnj - Muym] = 0. 



(35) 



However if we impose DBC's on the A5 components then, in the unitary gauge, ( [60| ) and ( [67] ) imply that 
one must impose DBC's on the pseudo-scalars, 0vKz|^_ 
(|59l) and SGE^, 



= 0. To see that this is consistent we substitute 



M^Vr^g5 (/i^^'^VS^^) ^ M-^]^a-^b-'d,(a'b-'Ml^^) 



A'pw.z) ,{") 



(36) 



where the + refers to tti^ and 7:2 respectively, into the boundary conditions of tt,- and note that the second 
term in ( [36] ) will vanish by the DBC's. While if we impose DBC's on the Higgs VEV, h, then the first term 
will always vanish. Alternatively if we impose non DBC's ([8]) then we must impose impose non DBC's on 
the Higgs components ( [35) . 



The upshot of all this is that if NBC's are imposed on the A^ field then the pseudo-scalars must have DBC's, 
regardless of the boundary conditions of the Higgs. This should be intuitive if one considers the zero mode 
of the field which will be a massless 4D gauge field and hence will have 2 transverse degrees of freedom 
after compactification. By the same arguments, that we began this section with, then clearly the degrees of 
freedom would not match if the physical pseudo-scalars also gained a zero mode. 

'To be explicit, the gauge field transforms over the orbifold as A^{r) — > A^(-r) = PA^{r)P'' and A^{r) — > A^i-r) = -PAs(r)P'' . 
Where P is a unitary matrix and the relative minus sign, which preserves gauge invariance under Am — » Am - jc^m©, will ensure 
opposite BC's for A,, and Aj. 
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5.2 The Pseudo- Scalar Masses. 



For gauge fields propagating in the pseudo-scalar masses and profiles are found by solving (57i and 



dsfn 



i<Pw,z) 



M' 

'Mw,z 



m; 



/2 



+ 1 _4A" + 8A'2 + - 2^ - 4A'^;;^ - M^,^ + e^Vi^^-^'^ 



where ' denotes d^. Typically this equation cannot be solve analytically and hence we must work with 
numerics. The first pseudo-scalar masses have been plotted in figure |5] For the RS model (|4]) and Q, this 
gives 



M' 



w,z 



„akr 



M 



2k + ak- 



Be 



-akr 



W,Z 



,akr 



+ Be 



-akr 



and 



M 



w,z 



M 



w,z 



{A + a)k^ + Aak^^ 



akr 



Be 



-akr 



akr 



+ Be 



-akr 



and so as a —> then ( 37 ) will reduce to the same equations of motion of a gauge field of mass Mwz ( 54 1, i.e 
the mass eigenfunctions will be approximately root one Bessel functions. Also as either a ^ 00 (i.e. a brane 

M' 

localised Higgs) then 00 and the pseudo-scalars will become infinitely heavy. Likewise if g ^ 

(i.e. the gauge fields become decoupled from the Higgs) then cp^f^z and the field should be considered 
unphysical. It is also clear from figure [5] as well as the above equations, that as a is increased the mass of 
the pseudo-scalars will always increase. Hence, at tree level, the pseudo-scalars will always be heavier than 
the first KK gauge field. This feature further adds to the ability of the LHC to rule out this model. 



5.3 The Pseudo-Scalar Couplings to Fermions. 



We can now move on to look at how the pseudo-scalars couple to fermions. Again in the interests of 



generality we shall work with the generic metric (43 1. The Higgs sector and gauge sector couple to the SM 



fermions via the Yukawa couplings and the fermion kinetic terms 

- - b - - - - \ 

-Ty''D„*F + /^y^Ds^ + i-Xy^D^X + iXy^D^X - bYo'^^X - bYueab'^''^^''X + h.c. , (38) 



where ^ 3 |L, Q} and X 3 {e, u, d), as in section [3] The first and third y^D^ terms are just the usual gauge 
fermion interactions. The second and fourth y^Ds terms will give rise to pseudo-scalar currents interacting 

with the Aj components, for example. With - | q ^- j then 

where A^ are the fifth components of the SU(2) gauge bosons, as in the appendix. If we remember that, after 
compactification over the orbifold, only the left handed components of ^ and the right handed components 
of X will gain a zero mode then clearly the A^ component will only couple to currents containing one SM 
fermion and one KK fermion or alternatively two KK fermions. Hence at tree level, the pseudo-scalars will 
only couple to SM fermion currents via the Yukawa couplings to the Higgs components, ;r,. Bearing in mind 
that the larger the value of a, the more the pseudo-scalars will be dominated by the ni components and hence 
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Figure 5: The masses of the first Z pseudo scalar for the RS model (black), as well as the modified metrics 
with V = 10 (blue), v = 5 (green) and v = 3 (red) and A:A = 0.5 (dash-dot line), kA = I (solid line) and 
kA = 1.5 (dash-dash line). 



we can see our first indication that pseudo-scalars will have more significance for SM processes, for larger 
values of a. 

At next to leading order, the phenomenology becomes more complex and a proper investigation is beyond 
the scope of this paper. Although it is worth commenting that the Higgs particle will couple to the pseudo- 
scalars via the l^l"* term and so one would anticipate that the Higgs mass will be sensitive to the pseudo-scalar 
masses. Hence if a is very large and the pseudo-scalars are very heavy, then one risks reintroducing a little 
hierarchy problem. 



Moving onto look at the Yukawa couplings and looking just at the quarks, 
S 

J V2V 

-Yu U(^){h + H - ini)u(^) - Yu d(^)T:-U(j^) + h.c. 



= J d^x - Yd U(^)n+d(j^) - Yd d(if,){h + H + ini)d{x) 



(39) 



where n+ = 4=(7ri ± 1712), we have also split Q = i "^^^ | and the {ip,x) subscripts refer to whether the 
fermion originated from a SU(2) doublet or a singlet. To find the effective coupling to the SM fermions, we 



can now plug in ( [59] ), ( [66] ), carry out a partial integration and use the DBC's of (pw,z to get 

) 



Xeff. C 

V2 



m„ 



(<l>w,z) 2 



. Y^^'^^SJ^'^l (40) 
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For the spaces considered in this paper this gives 

where N^^^^ are the fermion normalisation constants. At first sight, one may be alarmed to see that this 
effective Yukawa coupling appears to be dependent on the gauge coupling, g. However this is not the case 



since the normalisation constant of f^^\ obtained from (61 1, will contain a factor of g'^. Likewise the 
effective coupling will be independent of Nh at tree level. This is an interesting coupling for a number of 
reasons. Firstly it gives a good example of a scenario in which relatively small changes in the geometry can 
result in significant changes in the phenomenology. In the RS model, the 'e^'"' factor' in the h' term will 
exactly cancel with the 2A'h term and the effective coupling will be given by 



r 



YgNhk' j (1 - 2cl)(1 + 2cr) 1 



xjdr /i^"'^ (^e^~M)kr + CRXe"'"' + Be-"'"') - aie"'"' - Be'"'"'))) . (42) 

This coupling is very sensitive to the value of a, as can be seen in figure |6] For large values of a the pseudo- 
scalars would become strongly coupled and one would loose perturbative control of the theory. On the other 
hand, this cancellation does not occur in spaces with a modified metric. The resulting effective coupling has 
additional terms which result in the coupling being significantly less a dependent. 

The second interesting feature of this coupling is that there exists special fermion positions (when lA'h - 
h' - (cl - CR)kh K 0) for which the pseudo-scalars tree level coupling, to the SM fermions, is significantly 
reduced. This is analogous to the case of the KK gauge fields coupling to SM fermions, when cl - -cr - 
0.5. These fermion positions typically occur when the fermions are heavily peaked towards the IR and for 
many of the modified metrics considered, this occurs for bulk mass parameters not plotted in figure |6] 

Clearly there is significant room for further work concerning the phenomenological implications of these 
pseudo-scalars. However by going further in this direction we risk digressing too far from the central 
focus of this paper, that of investigating the fermion mass hierarchy and looking for possible motivations for 
considering a small value of a to be more plausible than a large value. Nevertheless, since we have suggested 
that the observation, or lack of observation, of such scalars could play an important role in excluding this 
model it is worth making some comments concerning their phenomenology. In addition to the Yukawa 
couplings, the pseudo scalars will also couple to the y, W and Z gauge fields, the KK gravitons, the radion 
and the Higgs itself. Also, due to their large masses, one would anticipate that they could decay via a KK 
particle. Hence any study of the production and decay of such particles would be quite involved and should 
really be conducted as a separate piece of work. 



Returning to the focus of this paper, after one has diagonalised the fermion mass matrix ( 19 1, this model 
would give rise to the possibility of tree level pseudo-scalar mediated FCNC's. The constraints from these 
FCNC's have the potential to force the KK scale to a level at which it becomes questionable whether one 
has resolved the hierarchy problem. The extent to which such FCNC's would be suppressed is determined 
by the size of the effective Yukawa couplings as well as the size of its misalignment with the fermion masses 
or equivalently its universality with respect to fermion position. This has been plotted in figure [6] It is 
found that, when the fermions are peaked towards the UV (c > 0.5), although this coupling is significantly 
misaligned and non-universal, it is also very small and hence such FCNC's would be heavily suppressed. 
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(a) The effective coupling of the first 0z field to the SM fermions assuming 7=1. 




(b) The effective Yukawa coupling relative to the fermion mass. The less universal this value 
is, then the greater the pseudo-scalar mediated FCNC's will be. 

Figure 6: The couplings of the Z pseudo-scalar to SM fermions, with bulk mass terms c, for the RS model 
(black) and modified metrics with v = 10 (blue), v = 5 (green) and v = 3 (red). Here we fix = 1 
and consider a = 0.01 (solid line), a = 0.5 (dash-dash line) and a = 1.01 (dash-dot line). D. = 10^^ and 
Mkk = 2 TeV. 

However, for the RS model, as a is increased the size of this coupling increases and it rapidly becomes 
potentially problematic. Surprisingly, this effect is not as severe in spaces with a modified metric, for 
reasons discussed above. 
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6 Implications for Flavour Physics. 



Of course, this model would not just receive constraints from pseudo-scalar mediated FCNC's. It is well 
known that models with warped extra dimensions, that seek to explain the fermion mass hierarchy, suffer 
from severe constraints from flavour physics (see for example |[34ll42ll36ll43ll44ll45l ). These constraints 
fall into three broad categories: 

• Those that are dominated by the modification of a SM coupling. These include, for example, correc- 
tions to the Zbb vertex as well as rare leptonic decays such as // ^ eee. 

• Those that are dominated by the tree level exchange of a KK particle. For example, a particularly 
stringent constraint arises from ex that receives large contributions from the exchange of a KK gluon. 

• Those that arise at next to leading order via additional penguin diagrams. In particular, constraints 
from /J. ^ ey ||46l|47l and b ^ sy ||48]| transitions. 
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Figure 7: Coupling of the Z zero mode to fermions at positions c - ci = -cr for the RS model (black) 
and the IR modified spaces with = 1 and v = 10 (blue), v = 5 (green), v = 3 (red). We use the Higgs 
exponent of a = 0.01 (solid lines), a = 1.01 (dash-dash line) and a = 5 (dash-dot line). In the interests of 
making a fair comparison we take a KK scale such that the mass of the first KK gauge field is the same, i.e. 
for RS Mkk = 2 TeV, for v - 10 Mkk = 2.15 TeV, for v - 5 Mrk = 2.70 TeV and for v - 3 Mrk = 4.59 
TeV. Also Q. - 10^^ 

While it certainly lies beyond the scope of this paper to conduct a full and thorough investigation of flavour 
in models with a bulk Higgs, we are in a position to see how the first two categories of constraints would be 
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effected. It is more difficult to estimate the effect on tlie next to leading order processes due to the number 
of conflicting factors contributing to these constraints. In particular there would be modifications to both the 
KK masses and the couplings plus additional diagrams coming from the pseudo-scalars and the KK Higgs 
bosons. 



6.1 Modifications to the SM Fermion Couplings. 



Of particular importance, to many constraints from flavour physics, is the modification to the Z fermion 
coupling, which has been plotted in figure|7] This coupling not only adds a direct constraint on the KK scale, 
via rare lepton decays and the partial decay width of Z ^ bb, it also constrains the extents to which a split 
fermion scenario is allowed. Such a modification occurs because, after spontaneous symmetry breaking, 
the Z zero mode is not flat but modified in the IR. This gives rise to non-universal couplings to fermions 
in different locations. As one would expect, the flatter the Higgs VEV the smaller this modification will 
be. Hence the smaller the value of a the smaller these constraints will be. It is this effect which is partly 
responsible for the reduction in the constraints from rare lepton decays seen in ||2T1|491 . Although it is also 
found that modifications to the geometry typically give an enhancement of the deformation of the Z zero 
mode, relative to AdSj. 



The fermion zero mode profiles will also receive corrections arising from the mixings between the SU(2) 
singlets and doublets, see (21 1. These corrections can also give a sizeable correction to the Z ^ bb vertex 
pij. The size of these corrections are determined by the size of F^*]^"^ for « ^ 1. Whether this correction 
increases or decreases, as one increases a, is not entirely clear since it is dependent on whether the fermion 
zero mode profiles are peaked towards the UV brane or towards the IR brane. Hence it is sensitive to the 
preferred fermion positions (see the following section). It is also sensitive to how quickly the Yf^^^ decreases 
as one increases n. It is straight forward to check that Yf^^ will drop away more quickly for smaller values 
of a. In other words, as one shifts the Higgs VEV away from the IR brane one reduces the mixing between 
fermion zero modes and the higher KK fermion modes. 



Another source of flavour violation, in these models, would be the Higgs particle couplings to the SM 
fermions. Since the Higgs particle profile (16 1 is not the same as the Higgs VEV (j9]l, after diagonalising 

The size of 



( 19 1, the resulting Higgs particle effective Yukawa couplings will not be flavour diagonal II50II 
this effect will be determined by the difference between the Higgs VEV and the Higgs profile which, after 

expanding ( 16 1, is found to be ~ O ( 



2a 



Hence this is the one counter example considered in 
which reducing a would increase the size of the correction. However in practice it is found, for all spaces 
considered, that when mj, «; Mkk this effect is quite small. 



6.2 FCNC's from the Exchange of a KK Particle. 



As already mentioned, models with warped extra dimensions and bulk fermions suffer from stringent con- 
straints from FCNC's arising at tree level via the exchange of a KK particle. In particular from - 
mixing and the observable e/f. Such FCNC's occur at tree level because the KK gluon has a non-universal 
coupling to fermions with different bulk mass parameters. Models with warped extra dimensions and a 
large warp factor have a natural suppression of such FCNC's referred to as the RS-GIM mechanism (see 
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Figure 8: Plotted is the relative coupling of a KK gluon for the RS model (black) and the modified metrics 
with V = 10 (blue), v = 5 (green) and v = 3 (red). As well as with /cA = 1.5 (dash-dash line), kts. - \ (solid 
line) and ^A = 0.5 (dash-dot line). These couplings are independent of a at tree level. The relative difference 
in the couplings between fermions peaked towards the UV and a fermion with a bulk mass parameter of 
c = 0.6 has been plotted on the right. Q = 10^^ and Mrk = 2 TeV. 



for example l'34'l). In particular fermions peaked towards the UV (c > 0.5) have approximately universal 
couplings to the KK gluon, see figure [8] The problem really arises because typically it is found that the 
heavier quarks have to sit quite far towards the IR and hence have non-universal couplings. 



The effect on this constraint, of allowing the Higgs to propagate in the bulk, has been studied for the RS 
model pO'l, spaces with a modified metric 1 17 1 and spaces with a soft wall [21 J. All three papers found a 
sizeable reduction in the constraint relative to that of the RS model with a brane localised Higgs. It is easily 
checked that the couplings of a KK gluon propagating, in the modified metrics, are not significantly modified 
relative to that of the RS model, see figure [8] This coupling is also found to be independent of a. Hence it 
is suspected that the primary reason for this reduction in the 6k constraint is due to a shift in the preferred 
fermion positions. In particular the heavy quarks can sit further towards the UV, where the coupling is more 
universal. When referring to the preferred fermion positions we are referring to the bulk mass parameters, 
cl,r, that give the correct masses and mixing angles assuming anarchic order one 5D Yukawa couplings, Y. 



However none of the above papers investigated the a dependence of this shift in the fermion position. A 
priori it is not obvious how the fermion positions are effected since it is based on two conflicting factors. 



Notably the 'maximum fermion mass' decreases as a is reduced (see figure 3(b) I but the gradient of the 



slope in figure 2(b) also reduces. Hence this must be calculated explicitly. To do so we minimise a function 
that inputs the nine quark bulk mass parameters (cl, and c^) and computes the value for the mean 
quark masses, mean CKM mixing angles and mean Jarlskog invariant, taken over 5000 anarchic Yukawa 
matrices with ^ < |F| < 3. The quark masses are run down to the mass of the first KK gluon from the 2 GeV 
values on. 



Mu ^ 2.4 ± 0.7 MeV 
md = 4.9 ± 0.8 MeV 



Mr 



1.29!°;?] GeV 



m, = 100!^° MeV 



rrit = 173 ± 1.5 GeV 
mh = 4.2!°;!^ GeV 
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and mixing angles used are IST] |5T]| 

y„,s = 0.2254 + 0.00065 V,/, - 0.0408 ± 0.00045 - 0.00376 ± 0.0002 7 = 2.9l!° x 10"^ 

Due to the size of tlie parameter space, such a minimisation routine will typically find a local minimum 
and not a global minimum. Hence this is repeated 200 times with random initial guesses in the ranges 
CL = [0.66 ± 0.1, 0.58 ± 0.1, 0.4 ± 0.2], 4 = [-0.6 ± 0.1, -0.51 ± 0.1, 0.2 ± 0.5] and 4 = [-0-66 ± 
0.1, -0.64 + 0.1, -0.59 ± 0.1]. Of these 200 configurations, the 40 best^^ values are taken and the mean 
bulk mass parameters have been plotted in figure|9]for difl'erent values of a. Although some prior knowledge 
of warped extra dimensions has been used in choosing the initial guesses, it is still found that all the bulk 
mass parameters do still converge to a preferred value. 

Inevitably, due to the size of the parameter space, there are large numerical uncertainties associated with 
these preferred fermion positions. Also one could argue that a better analysis would include more observ- 
ables such as ex and the Z partial decay width, /?/,, in the fit. These would slightly shift the preferred values 
but it is suspected that the basic result would still holcj^ The result being that, despite these uncertainties, 
there is a generic trend for the fermions to sit further towards the UV for small values of a. This would 
lead to their couplings being more universal and hence result in a reduction in the constraints from not just 
6/f but a significant number of the AF - 2 and AF = 1 processes. One can see from the analysis made in 
ifnl I20I I2TI that relatively small changes in the fermion positions can result in quite significant reductions 
in these constraints. 

Finally it should be pointed out that another source of flavour violation in these models arises from Higgs 
mediated FCNC's arising from the exchange of a KK fermion. This has not been considered here since it 
was studied in ||50l and found to be largely independent of a. 



7 Conclusions. 



The primary focus of this paper has been to investigate the phenomenological changes, in RS type scenarios, 
as one changes the exponent of a bulk Higgs VEV. The motivation for doing so is that, when the exponent 
is close to two and one assumes order one Yukawa couplings, the range of fermion zero mode masses is 
in remarkable agreement with the observed fermion mass hierarchy. After first introducing the model, it 
was demonstrated that, when brane localised potentials exist, then non tuned Higgs VEVs will not be flat, 
but will be peaked towards the JR. Hence bulk Higgs scenarios still offer a potential resolution to the gauge 
hierarchy problem. However, models with a bulk Higgs, gain an additional positive contribution to the lOp 
term in the effective potential which may result in EW symmetry not being broken for large values of the 
Higgs exponent, a. Next the fermion mass hierarchy was considered and it was found that (assuming order 
one Yukawa couplings) the Dirac mass term of the fermion zero modes stretch from the EW scale to a factor 
of about " below it, see figure [s] for the calculated values. We then proceeded to demonstrate that a 
considerable number of the phenomenological constraints will be at a minimum for smaller values of a. In 
particular it is shown, for generic Higgs potentials and generic geometries, that the flatter the Higgs VEV 
the smaller the contribution to the T parameter. It is also shown that a large class of constraints from flavour 

'^The a dependence of the fermion position is dominated by tfie gradient of tfie slope and the position of the 'maximum fermion 
mass' in figure [2(b)l Including additional observables would not change this but would rather change the extent to which a split 
fermion scenario was favoured. For example the inclusion of would favour a shift in Cl towards the IR and towards the UV 
12H . While the inclusion of Rt would disfavour a split fermion scenario. 
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Figure 9: The mean and standard deviations for the 'preferred fermion' positions. When cl > 0.5 (< 0.5) 
the left handed fermions profiles will be peaked towards the UV (IR) branes and likewise when cr < -0.5 (> 
-0.5) the right handed fermions will be peaked towards the UV (IR). Here we have considered the RS model 
(black) and the modified metric with = 1 and v = 10 (blue) and v = 3 (red). The light quark positions do 
not demonstrate as significant an a dependence. Q= 10'^ and Mkk - 2 TeV. 



physics will be reduced for smaller values of a. This is related to three effects. Firstly for the RS model the 
size of the pseudo-scalar coupling to the SM fermion is enhanced as one increases a. This coupling is not 
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flavour diagonal and would give rise to pseudo-scalar mediated FCNC's. This effect is significantly reduced 
in spaces with a modified metric. Secondly the non-universality of the W and Z coupling to SM fermions is 
reduced for smaller value of a. Thirdly it is found that there is a general trend for the fermions to sit further 
towards the UV for smaller values of a. This would result in their couplings to KK gauge fields being more 
universal with respect to flavour. 

All of the above motivations seem to favour a small value of a and aie independent of the requirement of 
generating the correct fermion mass hierarchy. One also has a potential lower bound of a > coming 
from the Breitenlohner-Freedman bound, although this is not strictly applicable here. However the RS 
model is not a UV complete theory and it is plausible that attempts to embed such a model in a more 
fundamental theory could result in the Breitenlohner-Freedman bound being applicable. Another lower 
bound arises when one considers the 4D dual theory. The corresponding Higgs operator, O, would have 
a scaling dimension of 2 -i- a. In order to off'er a potential resolution to the gauge hierarchy problem one 
requires that the operator O^O is not relevant which in turn implies that a ^ [6]. Hence it is compelling to 
saturate these two bounds, for the largely phenomenological reasons given in this paper and claim that the 
optimal value of a is zero (or equivalently the optimal exponent of the Higgs VEV and scaling dimension 
of the Higgs operator is two). 

Even if one accepts that a should be close to zero, then more work is still required before the fermion mass 
hierarchy can be understood. In particular, models with a bulk Higgs offer no explanation of, for example, 
why the leptons are lighter than the quarks, or why there is a six orders of magnitude gap between the 
electron mass and the neutrino masses. We have also focused on the Dirac mass term. It was argued in (5], 
that in RS type scenarios constraints, from lepton flavour violation, favour Dirac neutrinos over Majorana 
neutrinos. However there is no symmetry forbidding a Majorana mass term. The validity of this possible 
explanation of fermion mass hierarchy is conditional on such mass terms being either forbidden or existing 
at the same mass scale as the Dirac term. There are a number of possibihties for achieving this but all would 
require an extension of this minimal model and so we do not consider them here. 

Further work is also required in order to fully understand the phenomenological implications of both modi- 
fying the geometry and the pseudo-scalars. With regards to the modification of the geometry, although here 
we have not fully explored the possible parameter space, we still find significant phenomenological changes 
relative to the RS model. In particular it is found that, while the couplings do not change significantly, 
the relationship, between the curvature and size of the extra-dimension and the KK mass eigenvalues, does 
change. It is believed that this is partly responsible for the reduction in the EW constraints lfT4l [TSl Il6l . 
With regards to the additional scalar degrees of freedom, assuming the LHC is able to reach the KK scale, 
then the prediction of a pseudo-scalar, for each W and Z', makes models with a bulk Higgs very falsifiable. 
Although again we must leave a full study of such scalars to future work. Nevertheless, it is believed that 
there are a number of reasons for considering models with a bulk Higgs, an interesting extension of the 
description of flavour that already exists in the RS model. 
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A Appendix 

In this appendix we shall derive many of the results used throughout this paper. Here we will work with a 
generic 5D warped space of the form, 

ds^ = airfrfdx^dxy - b{rfd?, (43) 

with r e [ruv, Hr]- Without loss of generality, b{r) can always be set to one with the coordinate transforma- 
tion r — > r = b(r)df. However, by not doing so, one can analytically express a greater range of spaces 
and also one can easily modify the following expressions in order to use a conformally flat metric. 



A.l Electroweak Symmetry Breaking with a Bulk Higgs. 



Let us begin by considering a bulk SU(2) x U(l) gauge symmetry in addition to a bulk Higgs, 3), 

5 = Jd'x VG|-^F^^Ff ^ - ^-BmnB'''' + |DmO|2 - y(0)J . (44) 

Where the covariant derivative is given by Dm - Om - ig^^'^" - i^g'^M, with the three SU(2) generators 
t" - (t"/2 and a = {1,2,3). After spontaneous symmetry breaking the Higgs would acquire a non zero 

VEV, (3)) = \ where h would be the solution of 

ds (a^b-^dsh) - a'^bd,s,Vm\<!,=h = 0. (45) 



V2 \ h(r) 



Clearly in order to actually break EW symmetry one must have a potential that does not admit /i = as a 
solution. This can be done by either using the bulk potential, V{(i>), or by imposing non-trivial boundary 
conditions. We can now expand 3) around h such that, 

V2 \ hir) + Hix, r) + im(x, r) ) ' ^^""^ 

and make the usual field redefinitions, 

^M-A^ {a'm + iAl) Zm = J— (gAl - s'Bm) Am = J— (g'A^ + sBm) . (47) 
V2' yfF+g^ iF+g^ 

The Higgs particle itself will satisfy 

a^bd^d^H - dsia^b'^dsH) + a'^bdq>Vm\q>=H = 0. (48) 
It is also useful to define the quantities 

Mz(r) . Vi^Mr ) M^(r).^. (49) 
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These, of course, should not be confused with the 4D W and Z masses. We must now include a gauge fixing 
term which is chosen in order to cancel the mixing between the 4D gauge fields, and the 4D scalars, 
and Ka- In particular, working in the gauge, we introduce the gauge fixing term 

-CcF. = {d^Z^' - ^b-' [d^ia^b-'Zi) + a^bMz7ri)f - ^ [d^A^; - ^b'' (dsia^b'^Al) - a^bMw7T2)f 

{dXi - (dsia'b-'A',) + a^bMwni)f . (50) 

If, for the moment, we just focus on the Z field then the equations of motion are found to be 



bUzy^ + ^d^id.Z')] + ^5 {a-b-'dsZ^) - a'bMlZ^ = 0, (51) 



d^d^Zs + a^Mzdim - a^Mz + a^MjZs - ^85 (b~^ [diia^b'^Zs) + a^bMzn^)) = 0, (52) 

df.ff'ni - a-%-^ds [a^b-^din^) - a'^b'^d^ [a^b'^MzZi) - a^b'^Mz (^J^s + a'^d^Vmh=„, 

+^Mzb~^ [dsia^b'^Zs) + a^bMzT:^) = 0. (53) 

To find the masses of the 4D Z bosons, one would expand Z^ in terms of orthogonal mass eigenstates, i.e. 
make a KK decomposition, Z^(a;, r) = fn^\r)Zj^\x) such that ^ drbfjp f^^^ - 6„,n- The masses would 
then be defined by the 4D equations of motion, 5^Zy^^ + |5^(5''Zy"^) = ni^^^^zjf\ and found by solving for 
the profiles, 

^5 (a^b- ' - a^bMy^^^ + bmf^ ^ - 0. (54) 

The W and photon profiles, fj,^^ and /^^^ are found using the same equation but with Mz replaced by Mw 
and respectively. 



A.2 The Gauge-Goldstone Boson Equivalence Theorem Cross Check. 



There are now two remaining scalar degrees of freedom. One will be eaten by the longitudinal polarisation 
states of the 4D massive Z fields and will be the unphysical Goldstone boson. The other will be a physical 
pseudo-scalar. To find the Goldstone bosons we simply note that such fields must have a completely gauge 
dependent mass. Such that in the unitary gauge (^ 00) the Goldstone bosons become infinitely heavy and 
hence should be considered unphysical. From ( [52] ) and ( [53] ) it is straightforward to see that the Goldstone 
boson should be given by 

0z = b-^ (85 (a^b-^Zs) + a^bMzKj) . (55) 

However in order for the gauge-Goldstone boson equivalence theorem to hold, in the 4D theory, one nec- 
essarily requires that the profiles of the Goldstone boson, f^^^\ should match those of the Z boson, 
Otherwise the 4D effective couplings would be different for the two fields and hence the scattering ampli- 
tudes would not be the same. This equivalence has been found in other extra dimensional scenarios, such 
as when a gauge field propagates in more than five dimensions ||52l . To check that the profiles do match 



a^bMz\53\ and noting that dsia'^b-^MjZs) - Mzdsia'^b-^MzZs) - a^b-^Ml{^)Zi = gives 



we must find the equations of motion for the Goldstone bosons. In particular, taking 85(0 b (|52|)), adding 

qZs) - Mzdsia^b-'MzZs) - a^b-'M^C^ 



d^d^'gz - ^^^^^'^ ds {a%-'d5h)n3 + a%Mzd^Vm^=,, - ^b'^ds (a^b-'d^Qz) - ^a^M^Qz - 0. 
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For all Higgs potentials, considered in this paper, the second and third gauge independent terms will cancel 



using (45 1. One would anticipate that such a cancellation would occur for all Lorentz and gauge invariant po- 
tentials, but this is not verified here. With such a cancellation, one can proceed to make a KK decomposition 
Qz = Y.nff^\r)QfiX) such that d^&'Qf - -rf'^'^^z^ and the profiles will be given by 



(56) 



in agreement with ( 54 1. 



A.3 The Physical Pseudo-Scalars. 



To find the remaining physical degree of freedom we must take a combination of (52) and (53 1 such that 



the mass term is gauge independent. This can be done by adding ( |52[ ) to dsiM^ {53 1) which, after a little 
algebra, gives 



We have again used the cancellation ds{a^b'^ d^h) + a^bd(SfV{Q>)\(^=ni = 0, as well as defining 



(57) 



(58) 



This equation agrees with equivalent expressions derived previously in |fT5ll4ni . In order to find the effective 
action for these pseudo-scalars, we again make a KK decomposition, (f>z{x, r) = f!f^\r)<Pz^{x), such that 



d^&'ip^^^ = -m,f^^ '^(p^^K This allows us to invert \55) and (|58| to get 



=2 



[ M~^a~-b~^8,(a^b~^Mlfl?^^)^^^^ Mz_g^^,^ 

(0z)2 



m 



'PZ T ^ 



ml, 



z)2 



m. 



i@z)2 ^Z 



(59) 



(60) 



For the bulk of this paper we shall work in the unitary gauge, in which the Goldstone bosons are infinitely 
heavy and can be neglected. The final step is to ensure the mass eigenstates are canonically normalised 
which is achieved with the orthogonality relation 



dr a'b 



Mi 



Mz) A</>z) _ 
Jn Jm 



(61) 



A.4 The W Boson. 



It is straight forward to repeat this analysis for the W field, where it is found that the two Goldstone bosons 
are given by 

§\^ = b-^ds{a%-''A^^-a%MwTi2 and §1^ = b-^ds[a%-^Afj + a%Mwni. (62) 
While the two physical pseudo-scalars are found to be 

'Pw = d5{M-y)n2)-A\ and cp^ = dsiM-^'m) + A]. (63) 
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In practice there will be a further mixing such that the charged pseudo-scalars are given by combinations of 
n- = -^(^1 ± 1^2) and = -^i^l + '^5)- After making analogous KK decompositions, as those used for 
the Z field, the equations of motion are analogously 



The effective action can be found by substituting for 



J ,2 



^1,2 



m^a D o^\a D m^Jn j 2.1(«) ^wL£_^2,\ (n) 



i<l>w)2 



m 



w 

W-fn A, 2 in) 



m 



i<t>w)2 



(&w) 



m. 



(&w)2 



m. 



^1,2 (n) 



(64) 
(65) 



(66) 
(67) 



A.5 The Fermion Profile. 



For completeness, we shall briefly derive the widely known expressions for the fermion profiles. Beginning 
with the action for a massive Dirac spinor, *F, in 5D 

S= J^j5^VG(;Tr^VM'F-Mvp4"F) (68) 

where the Dirac Matrices in curved space = E'^j^ are related to Dirac matrices in flat space, y^, by 
the fiinfbein - diag(a^\fl;~^a~^a~',Z7"'). While the covariant derivative, - Dm + ojm, includes 
the spin connection ojm - (^^"'^sa 757^,0). Splitting ^' = i^l + if/R such that iysifj^R = +^l,r and making 
the KK decomposition ifji^R = a-^f^f"-"\r)ilf^l^l{x) such that {iy^df, - m^f)'i'^"^ - 0. This then yields the 
coupled equations of motion 

+ bM^f!f^' = -mf/^^ and - ^5/^ + ^M^p/^^ ^ -mf /f (69) 
a a 

It is then straight forward to solve for the fermion zero mode profile 

^^iMff^^^. (70) 
^filriexp{±2£dfb{f)M^) 

It is common to parameterise the bulk mass term such that Mvy = -ck where c ~ 0{\). 
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